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1.0  INTRODUCTION 


This  report  presents  the  results  of  a  computer-based  wave  refraction 
analysis  which  was  executed  in  conjunction  with  the  design  of  the  four 
marine  structures  of  the  U.S.  Navy  East  Coast  Air  Combat  Maneuvering 
Range,  located  offshore  at  Kitty  Hawk,  North  Carolina.  Figure  1-1  shows 
the  general  layout  of  the  towers. 

The  structures  will  consist  of  a  three-pile  jacket  (template)  with 
equi laterally-spaced  legs  through  which  steel  piles  are  driven  into  the 
seabed.  Jackets  are  then  secured  to  the  piling  by  welding  shim  plates  in 
the  annulus  between  the  jacket  leg  and  piling  at  the  top  of  the  jacket 
legs.  Following  this,  a  superstructure  consisting  of  an  equipment  deck 
and  an  upper  deck  is  attached  to  the  piling  above  the  top  of  the  jacket. 

Orientation  of  the  structures  will  be  such  that  the  side  of  the  plat¬ 
form  with  solar  panels  will  face  due  south.  This  places  the  side  of  the 
structure  with  the  proposed  boat  landing  to  the  northeast.  Site  locations 
of  the  structures  are  presented  in  Figure  1-2  and  will  be  as  follows: 

Site  #  Latitude  (N)  Longitude  (W)  <“> 


1 

35°  56.99' 

75° 

16.00' 

81 

15.80 

2 

36°  13.60' 

75° 

14.99' 

93 

24.00 

3 

36°  03.87' 

74° 

59.00' 

105 

31 .70 

4 

35°  47.19' 

75° 

05.71' 

105 

20.00 

The  goal  of  the  wave  refraction  analysis  was  to  determine  if  the 
topographic  refraction  effect  on  deep  water  waves  propagating  toward 
shore  would  produce  anomalously  high  local  v^ave  conditions  at  the  pro¬ 
posed  tower  construction  sites. 
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2.0  ANALYSIS  CONSTRAINTS 


“  V  r  fc.  ••  wj 


The  wave  periods  considered  in  this  analysis  are  10.0,  11.5,  and 
13.0  seconds.  It  was  assumed  that  the  waves  likely  to  have  the  greatest 
impact  on  structure  survival  will  fall  into  this  range  of  periods.  Deep 
water  wave  directions  (9)  considered  for  each  wave  period  of  concern  are 
the  following  azimuths  (defined  as  forward  direction  of  wave  propagation 
measured  clockwise  from  north):  210°,  225°,  240°,  255°,  270°,  285°,  300°, 
315°,  and  330°.  Thus,  a  total  of  27  deep  water  wave  conditions  were  con¬ 
sidered  in  the  analysis. 


•  • 


•  • 


•  •  • 


3.0  DIGITIZED  DEPTH  FIELD  WITH  GRID  COORDINATES 


To  proceed  with  the  wave  refraction  analysis,  it  was  necessary  to 
produce  a  digitized  depth  field  with  grid  coordinates.  The  depth  grid 
for  this  analysis  was  based  on  the  bathymetric  chart  of  the  region 
prepared  by  V.  Goldsmith  and  C.  M.  Sutton  of  the  Department  of  Geological 
Oceanography,  Virginia  Institute  of  Marine  Science.  Included  in  this 
"gridded"  depth  field  was  the  area  bounded  by  latitude  35°  20.5'N  to 
37°  00.0'N  and  longitude  74°  35.0'W  to  75°  24.5‘W  (see  Figure  1-2),  The 
grid  interval  is  0,5  minute  in  angular  distance  (both  latitude  and 
longitude),  whereas  the  origin  of  the  grid  coordinates  (1,1)  corresponds 
to  35°  20.5 ‘N  and  75°  24.5 'W,  A  right-hand  coordinate  system,  positive 
toward  east  and  north,  was  established  which  uses  100  by  200  grid  points. 
Figure  3-1  presents  the  "gridded"  depth  field  which  is  keyed  to  the 
geographical  area  of  interest  and  site  locations  presented  in  Figure  1-2. 
Depth  information  was  extracted  from  bathymetric  chart  at  every  third 
grid  point.  Depth  data  for  the  remaining  grid  points  was  obtained  by 
employing  an  interpolation  scheme.  Thus,  a  smooth  grid  coordinate  system 
for  the  depth  field  of  the  entire  area  of  concern  was  obtained.  The  grid 
coordinates  of  construction  sites  #1,  #2,  #3  and  #4  were  (18,73.98), 
(20.20,  107.20),  (52.00,  87.74)  and  (38.58,  54.38),  or  approximately 
(18,74),  (20,107),  (52,88)  and  (39,54),  respectively.  (Note:  Logi- 
tudinal  grid  points  are  given  first.) 


Note  that  the  "gridded"  depths  near  or  at  the  sites  may  differ 
slightly  from  recently  collected  soundings.  This  is  due  to  the  combined 
effect  of  depth  interpolation  to  derive  a  smooth  bathymetric  field  and 
the  intrinsic  inaccuracies  of  a  bathymetric  chart  based  on  a  variety  of 
sounding  sources.  This  discrepancy  is  of  little  consequence  to  the 
accuracy  of  the  wave  refraction  analysis. 


4.0  COMPUTER  ANALYSIS 


The  numerical  wave  refraction  analysis  involved  the  execution  of 
three  major  computer  programs.  The  first  program  performs  the  calcula¬ 
tions  of  all  relevatjt  parameters  along  an  individual  ray  on  the  global 
coordinates  for  a  given  wave  period,  initial  direction,  and  starting 
point.  Following  execution  of  the  first  program,  a  second  computer  run 
performs  a  correction  of  the  refraction  coefficient  near  and  at  a  wave 
"caustic".  A  caustic  point  occurs  where  two  adjacent  rays,  originally  a 
differential  distance  apart,  intersect.  Consequently,  the  ray  separa¬ 
tion  factor  becomes  zero  at  that  point.  According  to  the  conventional 
geometric  optic  ray  theory,  which  is  employed  in  the  first  program,  the 
refraction  coefficient  (reciprocal  of  the  ray  separation  factor),  and 
therefore  the  wave  height,  will  become  infinitely  high  at  the  caustic 
point  and  unrealistically  high  in  the  region  of  caustic  influence.  Thus, 
a  second  computer  run,  to  correct  for  this  error,  has  been  executed. 

The  third  computer  program  is  used  to  present  numerically  the  char¬ 
acteristics  of  refracted  waves  in  terms  of  the  wave  height  amplification 
factor  and  the  directional  change  from  the  initial  wave  direction  for 
each  grid  point.  This  is  accomplished  by  employing  a  weighted  average 
scheme  according  to  the  distances  from  the  grid  point  of  interest  to  the 
surrounding  ray  points. 
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5.0  RESULTS 


The  results  of  the  computer  analysis  for  the  four  sites  of  concern 
are  given  in  Table  5-1.  There  are  nine  directions  indicated  with  three 
periods  for  each  direction;  thus  there  will  be  twenty-seven  cases  for 
each  site.  For  each  case,  an  amplification  factor  and  a  refracted  wave 
direction  have  been  given.  The  refracted  wave  direction  is  the  direction 
that  a  wave  travels  at  the  site  after  it  experiences  the  effect  of  bottom 
topography.  The  amplification  factor  is  based  on  a  unit  amplitude  wave 
in  deep  water  being  propagated  shoreward.  For  example,  if  the  amplifica¬ 
tion  factor  at  the  site  for  a  particular  period  is  0.7,  it  infers  that 
the  wave  height  is  70  percent  of  the  deep  water  height. 

Examining  the  refracted  wave  direction  data  for  the  four  sites  (Fig¬ 
ures  5-1  through  5-4),  it  can  be  seen  that  waves  moving  shoreward  from  the 
northeast  (propagating  southwest  210°  to  255°)  will  refract  in  a  clockwise 
direction.  Waves  moving  shoreward  from  the  southeast  (propagating  north¬ 
west  285°  to  330°)  will  refract  in  a  counterclockwise  direction.  .  Eastward 
propagating  (270°)  waves  show  the  smallest  refraction  effect.  These  re¬ 
sults  are  in  keeping  with  basic  wave  refraction  theory. 

The  amplification  factors  associated  with  the  four  sites  (Figures  5-5 
and  5-6)  show  some  variation.  Sites  #1  and  #2,  which  are  closer  to  shore, 
have  significant  amplification  for  a  number  of  directions  (see  Figure  5-5). 
For  site  #1,  amplification  of  deep  water  wave  height  takes  place  at: 

a.  11.5  Sec  period,  225°,  factor  =1.1 

b.  10.0  Sec  period,  255°,  factor  =1.2 

Site  #2,  Figure  5-5,  shows  amplification  at  a  11.5-second  period,  with 
all  directions  greater  than  290  degrees.  The  site  #2  amplification  fac¬ 
tors  have  great  significance  when  one  considers  the  hurricane  generated 
seas  that  attack  the  coast  from  the  southeast. 
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Figure  5-1.  Refraction  Correlation  Diagram  (Site  #1 ) 
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DIRECTION  OF  PROPAGATION  AT  SITE 


Figure  5-2.  Refraction  Correlation  Diagram  (Site  #2) 
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Figure  5-4.  Refraction  Correlation  Diagram  (Site  #4) 
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A  geographic  presentation  of  the  amplification  factors  and  refrac¬ 
tion  directions  is  given  in  Figures  5-7  through  5-31.  The  overall  matrix 
represents  a  5  x  5  section  of  the  model  grid  centered  on  the  grid  point 
which  contains  a  particular  site.  Each  grid  point  in  the  Y  and  X  direc¬ 
tions  is  separated  by  0.5  minute  of  latitude  and  longitude,  respectively. 
The  elements  of  the  5  by  5  geographic  matrix  are  then  broken  down  into 
3  by  3  matrices.  The  3  by  3  matrix  gives  the  amplification  or  refracted 
direction  data  as  a  function  of  the  initial  direction  of  propagation  in 
deep  water.  A  key  to  the  3  by  3  sub-elements  is  given  in  Figure  5-7. 

The  direction  information  is  in  degrees  where  negative  is  left  or  counter' 
clockwise;  thus,  a  -6  for  an  initial  direction  of  210°  would  give  a  di¬ 
rection  of  204°  at  the  grid  point  of  concern. 

It  is  expected  that  the  geographic  distribution  of  the  refraction 
data  will  be  useful  if  any  relocation  of  the  structures  is  contemplated. 
When  using  the  data,  one  should  assess  all  periods  —  10.0,  11.5,  and 
13.0  seconds  --  before  deciding  on  a  preferred  direction. 
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SITE  NO.  :J _  lATITUDE:  35°  56.99'N 

WAVE  PPPTnn-  13.0  Sec _  l  ONGITUDE:  75°  16.00'W 


Figure  5-8.  Amplitude  Coefficients  as  a  Function  of 
Direction  (Site  #1,  13.0  Sec) 
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Figure  5-9.  Directional  Deviation  from  Initial 
Propagation  (Site  -1,  13.0  Sec) 
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Figure  5-10.  Amplitude  Coefficients  as  a  Function 
of  Direction  (Site  #1,  11.5  Sec) 
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Figure  5-11.  Directional  Deviation  from  Initial 
Propagation  (Site  #1,  11.5  Sec) 
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Figure  5-12.  Amplitude  Coefficients  as  a  Function 
of  Direction  (Site  =1,  10.0  Sec) 
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Figure  5-13.  Directional  Deviation  from  Initial 
Propagation  (Site  -1,  10.0  Sec) 
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Figure  5-14.  Amplitude  Coefficients  as  a  Function 
of  Direction  (Site  **2,  13.0  Sec) 
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Figure  5-15.  Directional  Deviation  from  Initial 
Propagation  (Site  #2,  13.0  Sec) 
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Figure  5-16.  Amplitude  Coefficients  as  a  Function 
of  Direction  (Site  *2,  11.5  Sec) 
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Figure  5-17.  Directional  Deviation  from  Initial 
Propagation  (Site  ^2,  11.5  Sec) 
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Figure  5-18.  Directional  Deviation  from  Initial 
Propagation  (Site  42,  10.0  Sec) 
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Figure  5-19.  Directional  Deviation  from  Initial 
Propagation  (Site  #2,  10.0  Sec) 
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Figure  5-20.  Amplitude  Coefficients  as  a  Function 
of  Direction  (Site  #3,  13.0  Sec) 
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Figure  5-21.  Directional  Deviation  from  Initial 
Propagation  (Site  -3,  13.0  Sec) 
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Figure  5-22.  Amplitude  Coefficients  as  a  Function 
of  Direction  (Site  #3,  11.5  Sec) 
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Figure  5-23.  Directional  Deviation  from  Initial 
Propagation  (Site  i3,  11.5  Sec) 
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Figure  5-26.  Amplitude  Coefficients  as  a  Function 
of  Direction  (Site  #4,  13.0  Sec) 
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Figure  5-27.  Directional  Deviation  from  Initial 
Propagation  (Site  #4,  13.0  Sec) 
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Figure  5-28.  Amplitude  Coefficients  as  a  Function 
of  Direction  (Site  #4,  11.5  Sec) 
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Figure  5-29.  Directional  Deviation  from  Initial 
IM  Propagation  (Site  #4,  11.5  Sec) 
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Figure  5-30.  Directional  Deviation  from  Initial 
Propagation  (Site  #4,  10.0  Sec) 
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Figure  5-31.  Amplitude  Coefficients  as  a  Function 
of  Direction  (Site  #4,  10.0  Sec) 
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6.0  IMPACT  ON  ENGINEERING  DESIGN 


The  accepted  wave  design  criteria  for  the  four  sites  are  as  follows: 


Site  # 

50  Yr  Storm  Wave  Ht 
(Ft) 

Period 

(Sec) 

1 

60.3 

13.6 

2 

60.8 

13.6 

3 

61.3 

13.6 

4 

61.3 

13.6 

Extrapolations  were  made  of  the  wave  climate  data  to  derive  a  data  base 
for  fatigue  analysis.  It  is  noteworthy  that  wave  refraction  renders  these 
appraisals  more  conservative  in  the  case  of  sites  #3  and  #4.  This  is 
mainly  due  to  the  fact  that  the  wave  amplification  factor  for  sites  #3 
and  #4  remains  less  than  unity  for  all  periods  and  directions  considered. 

For  sites  #1  and  #2  the  picture  is  not  as  clear.  In  the  case  of 
site  #1,  amplification  is  expected  for  10.0-second  and  11.5-second  waves 
propagating  toward  the  site  from  the  northeast.  This  evidence  indicates 
that  the  fatigue  analysis  data  for  this  site  should  be  biased  somewhat 
upward  (probably  less  than  1%).  In  the  case  of  site  #2,  strong  ampli¬ 
fication  is  indicated  for  the  11.5-second  waves  coming  from  the  southeast. 
Particular  consideration  should  be  given  to  the  fatigue  impact  for  this 
case  because  amplification  occurs  for  all  southeast  directions  considered. 
(This  effect  is  probably  less  than  2%. ) 
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